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Abstract. For the first time the large-scale clustering and the mean abundance of galaxy clusters are analysed 
simultaneously to get precise constraints on the normalized cosmic matter density Q m and the linear theory 
RMS fluctuations in mass erg. A self-consistent likelihood analysis is described which combines, in a natural 
and optimal manner, a battery of sensitive cosmological tests where observational data are represented by the 
(Karhunen-Loeve) eigenvectors of the sample correlation matrix. This method breaks the degeneracy between 
Q, m and as- The cosmological tests are performed with the ROSAT ESO Flux-Limited X-ray (REFLEX) cluster 
sample. The computations assume cosmologically flat geometries and a non-evolving cluster population mainly 
over the redshift range < z < 0.3. The REFLEX sample gives the cosmological constraints and their la random 
errors of Q m — 0.341 + 



and as = 0.711 



H0.039 
-0.031 ■ 



Possible systematic errors are evaluated by estimating the 



effects of uncertainties in the value of the Hubble constant, the baryon density, the spectral slope of the initial 
scalar fluctuations, the mass/X-ray luminosity relation and its intrinsic scatter, the biasing scheme, and the cluster 
mass density profile. All these contributions sum up to total systematic errors of an m =+°"g^ and a ag =t.g 162- 
Key words, cosmology: cosmological parameters - X-rays: galaxies: clusters 



1. Introduction 

Several observational arguments suggest that we live in a 
geometrically flat Universe in a phase of accelerated cos- 
mic expansion (Riess et al. 1998, Perlmutter et al. 1999, 
Stompor et al. 2001, Netterfield et al. 2002, Pryke et al. 
2002, Sievers et al. 2002). Strong indications for a low 
cosmic matter density come from, e.g., the abundance of 
galaxy clusters (e.g., Carlberg et al. 1996, Viana & Liddle 
1996, Bahcall & Fan 1998, Eke et al. 1998, Henry 2000, 
Borgani et al. 2001, Reiprich & Bohringer 2002, hereafter 
RB02), and from the large-scale distribution of clusters 
(e.g. Collins et al. 2000, Schuecker et al. 2001, hereafter 
Paper I) and galaxies (e.g., Percival et al. 2001, Szalay et 
al. 2001). A recent overview including also new results on 
gravitational dynamics, weak lensing and the baryon mass 
fraction in clusters of galaxies can be found in Peebles & 
Ratra (2002). 

With the present investigation we are improving the 
constraints on the matter density f2 m and the linear the- 



ory RMS matter fluctuations as in comoving spheres with 
8 hT x Mpc radius. We base our approach on observational 
estimates of spatial clustering and abundance of galaxy 
clusters. Simultaneous constraints on f2 m and as are de- 
rived under the assumption of a flat geometry of the 
Universe. Cluster abundances although almost insensitive 
to the cosmological constant SIa are known to be quite 
sensitive to fl m and actually one of the best ways of mea- 
suring as (White et al. 1993, Eke et al. 1996, Kitayama & 
Suto 1997, Borgani et al. 1997, Mathiesen & Evrard 1998). 
Our method uses the statistical information of the present 
epoch large-scale structure to characterize the cosmologi- 
cal model. In comparison to other methods of the assess- 
ment of the large-scale structure e.g. through the galaxy 
distribution, the analysis of the galaxy cluster distribution 
has the advantage that it relies on astrophysics mostly 
governed by gravitation, so that cluster masses and bias- 
ing factors can be obtained from first principles, especially 
when the clusters are selected in X-rays (e.g., Borgani & 
Guzzo 2001). 
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We use the ROSAT ESO Flux-Limited X-ray 
(REFLEX) galaxy cluster sample. In the REFLEX sur- 
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vey special care was taken to get a well-controlled sample 
selection and a high completeness (Bohringer et al. 2001). 
The sample consists of the 452 X-ray brightest southern 
clusters with redshifts mainly below z = 0.3, selected in X- 
rays from the ROSAT All-Sky Survey (RASS, Voges et al. 
1999) and is confirmed by extensive optical follow-up ob- 
servations within a large ESO Key Programme (Bohringer 
et al. 1998, Guzzo et al. 1999). The sample has been used 
to measure with unrivalled accuracy the cluster X-ray 
luminosity function (Bohringer et al. 2002), the spatial 
cluster-cluster correlation function (Collins et al. 2000), 
and its power spectrum (Paper I). 

In order to get accurate constraints on f2 m and a% 
we use several independent cosmological tests based on 
the dependence of the cosmic matter power spectrum and 
the volume on the values of the cosmological parameters. 
On the one side we measure the spatial fluctuations of 
galaxy clusters. This test is known to be especially sensi- 
tive to the shape of the matter power spectrum on large 
scales. On the other side we measure the average abun- 
dance of clusters. This test is known to be especially sen- 
sitive to the normalization of the matter power spectrum 
on small scales. Finally and independent from the mat- 
ter power spectrum, we count clusters to measure their 
abundance as a function of redshift and thus to measure 
the volume which again strongly depends on cosmology 
(redshift- volume test, e.g., Robertson & Noonan 1968). 

A natural way to combine the different tests is of- 
fered by the Karhunen-Loeve (KL) eigenvector analysis. 
The method was first applied by Bond (1995) to analyse 
cosmic microwave background (CMB) temperature maps 
and translated by Vogeley & Szalay (1996) to the case 
of the spatial analysis of galaxies. The KL eigenvectors 
\P are constructed in a way to obey orthogonality prop- 
erties which allow unbiased studies of fluctuations up to 
the largest scales covered by a survey. The method avoids 
artifical correlations introduced by the survey window be- 
tween different fluctuation modes, affecting under realistic 
conditions all power spectral analyses based on plane- wave 
expansions on large scales. Applications of the KL method 
to galaxy surveys can be found in Matsubara, Szalay & 
Landy (2000) and Szalay et al. (2001). Their KL tests are, 
however, still restricted to the analysis of the clustering 
properties of galaxies. 

To determine a precise eigenvector base VP which does 
not bias the final results, a more refined fiducial cosmo- 
logical model has to be assumed which is already close 
to our best model fits (the construction of the fiducial 
model is described in Sect.^J): a pressure-less spatially 
flat Friedmann-Lemaitre model, the cosmic matter density 
£lm = 0.3, the cosmological constant Ha — 0.7, the Hubble 
constant in the form h = -ffo/lOOkms -1 Mpc -1 = 0.7, the 
linear RMS normalization ag — 0.75, the spectral index 
of initial scalar fluctuations ng = 1.0, the baryon den- 
sity n b h 2 = 0.020, and the CMB temperature T C mb = 
2.728K. 

In our first KL study of the REFLEX sample we fol- 
lowed the approach of A. Szalay and collaborators, apply- 



ing the KL method to the analysis of three-dimensional 
spatial fluctuations (Schuecker et al. 2002, hereafter 
Paper II). The resulting constraints on the matter den- 
sity, 0.03 < 0, m h 2 < 0.19 (95% confidence interval), were 
found to be quite robust against major changes in the 
model assumptions. 

The present investigation completes our previous KL 
study by utilizing now both the spatial fluctuations of 
galaxy clusters and their mean abundance to get con- 
straints on £l m and erg. Some general features of the 
KL method relevant for the present work are outlined 
in Sect.|^. In Sect.|^ we briefly describe a model for the 
matter and cluster power spectrum already developed in 
Paper II. We further introduce a new theoretical model 
for the average cluster abundance which replaces the 
more phenomenological prescription used in Paper II. The 
REFLEX sample and the values of important model pa- 
rameters are discussed in Sect.|]. Sect.^ summarizes the 
cosmological constraints obtained with the REFLEX sam- 
ple. In Sect.^ we discuss the results and draw our basic 
conclusions. In the Appendix we derive equations used to 
relate cluster masses defined for different density contrast 
and background levels. 

2. General method 

Clusters of galaxies are counted in cells. The results are 
summarized in a vector D with elements containing the 
counts obtained within the individual cells. Similarily, a 
second vector N is introduced which contains the aver- 
age (expected) model counts. The latter are the diag- 
onal elements of a noise matrix N. After the equaliza- 
tion of the different noise in the cells by N' 1 / 2 (whiten- 
ing) and the KL transformation into the * vectorbase, 
the values of the cosmological parameters are estimated 
by minimizing the differences AB = B — < B > be- 
tween the KL coefficients B — ty T N~ x / 2 D determined 
by the observed sample, and the expected KL coefficients 
< B > = ^ T N +1 / 2 determined by a cosmological model 
(see Paper II), where T denotes the transpose of a matrix. 

More specifically, the columns of the unitary matrix "J 
are the KL eigenvectors of the so-called whitened corre- 
lation matrix R' = A -1 / 2 RN^ 1 / 2 , with the elements 
of R given by R tj = NiN^ + 5 %j Ni and £y = 
Wv~ Iv d^ri J v d 3 rj£(\ri — rj\). Here, £ is the cluster 
correlation function, and Vi the volume of the i-th count 
cell centered at the comoving coordinate T%. 

The elements of AB show random cell-to-cell fluctua- 
tions caused by both the Poisson count noise and the large- 
scale distribution of the clusters. For the correct weighting 
of AB these fluctuations are taken into account. Note that 
the analysis of the fluctuations alone already gives useful 
cosmological constraints as described in our first KL anal- 
ysis (Paper II). 

A simple way to perform the minimization which also 
utilizes the cosmological information contained in the ran- 
dom fluctuations between the cells is suggested by the ob- 
served frequency distribution of the KL coefficients. For 
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the REFLEX cluster survey we found (Paper II) that for 
large cell sizes, B follows a Gaussian distribution with 
high statistical significance (new results are given in Fig. [| 
of the present paper). A multivariate sample likelihood 
function of an assumed cosmological model should thus 
have the form 



C(B\x) 



exp [-§ trace (C _1 Z)] 
v /(2 7 r) M |detC7[ 



(1) 



with the sample covariance matrix Z = AB(AB) T , the 
model covariance matrix of the KL coefficients C = 
ty T R' modcl ty, and the parameter vector x specifying the 
cosmological model. We assume as usual that the maxi- 
mum of (]l|) determines the values of the model parameters 
which yield the highest probability of obtaining cluster 
abundances and fluctuations as large as observed. Note 
that in Paper II, Z was fixed by a phenomcnological back- 
ground model (independent from any cosmological param- 
eter) and only C was regarded as a free model-dependent 
variable. In the present investigation both Z and C are 
computed with observed and modelled data in a consis- 
tent manner (see Sec.||). 



3. Model specifications 

We briefly outline the model for the observed cluster power 
spectrum P b s (k) described in more detail in Paper II and 
a new model for the average cluster abundance N which 
replaces the empirical approach used in Paper II. 

The Gaussianity of the KL coefficients B measured 
with the REFLEX clusters on large scales together with 
the linearity of the KL transform suggest that the large- 
scale cosmic matter distribution is Gaussian as well and 
can thus be characterized by a (linear theory) matter 
power spectrum P(k). We assume adiabatic scalar pertur- 
bations and a cold dark matter (CDM) plus baryon fluid 
as suggested by the observed power spectrum of CMB 
anisotropics (see recent CMB measurements mentioned in 
Sect.0). The resulting P(k) is determined by the trans- 
fer function T x (k), the spectral index ns of the initial 
scalar fluctuations, and the normalization parameter <jg. 
The corresponding transfer functions of Eisenstein & Hu 
(1998) are used, providing a more accurate description of 
the linear theory power spectrum than fitting equations 
of the kind given in Bardeen et al. (1986). Note that the 
us as introduced here reflects the amplitude of the power 
spectrum without any non-linear corrections. 

Generally, P bs{k) is an average over evolving matter 
power spectra and clusters with different values of the 
biasing parameter b. We summarize this using the pre- 
scription of Matarrese et al. (1997) and Moscardini et al. 
(2000). Note that the prescription implies that the ob- 
served mass and redshift-averaged cluster power spectrum 
has the same shape, i.e., functional form as the underlying 
matter power spectrum. In order to illustrate the sensitiv- 
ity of the results to the biasing model we use for b alter- 
natively the model of Sheth & Tormen (1999) and Kaiser 



(1984). We verified the model of P bs(k) (see Paper II) 
with a large set of cluster samples extracted from the 
Hubble Volume Simulations (e.g., Evrard et al. 2002 and 
the references given therein). 

For the average cluster abundances the following model 
is used. For an unclustered distribution, the average num- 
ber of clusters expected in a small cell centered at Right 
Ascension a, Declination 5, and redshift z, per unit red- 
shift and solid angle u) is 



dN(a,S,z) cr 2 {z) 
dz dui H(z) 



dM 



dn(M, z) 
dM 



(2) 



with the usual notation H{z) = Hq E(z) and r(z) = 
c/H J Q E~ 1 (z')dz' . The average number of clusters ex- 
pected in a KL cell thus is Ni — J v dN. Recent CMB 
measurements strongly suggest a flat space (see Sect.|IJ) 
so that E 2 (z) = S! m (l + z) 3 + Q\. The mass limit M lim is 
obtained from the X-ray luminosity limit -L m j n at a given 
(a, S, z) and the empirical mass/X-ray luminosity relation 
of RB02 including its intrinsic scatter (see Sect. Q). They 
measured the masses within a radius with the average den- 
sity contrast 200 related to the zero-redshift critical den- 
sity p c , assuming an Einstein-de Sitter (EdS) geometry. 

For cosmological tests a mass function dn/dM is 
needed which is form-invariant under different cosmolo- 
gies. Jenkins et al. (2001) found from numerical simula- 
tions that 

dn(M, z) p(z) 1 

v ' = -0.315 Fy ' 



dM 



da(M) 



dM 



M a(M) 

exp{-|0.61-ln[cr(M,z)]| 3 - 8 } (3) 



obeys this invariance property when the masses M are 
defined by a spherical overdensity of 187T 2 related to the 
background matter density p — fl rn p c . The growth of cos- 
mic structures introduces a redshift-dependent a in the 
mass function, giving the standard deviation of the mat- 
ter density field smoothed by a spherical top-hat filter with 
the radius R = [2GM /^iJ 2 ] 1 / 3 . 

In order to combine the empirical mass/X-ray lumi- 
nosity relation and the mass function, two transforma- 
tions of cluster masses are necessary which are described 
in Sect.|i| and Appendix [A| Any change of the values of 
cosmological parameters thus changes both the observed 
mass/X-ray luminosity relation and the theoretical mass 
function. Note that the cosmology dependence of the ob- 
served mass/X-ray luminosity relation of RB02 results 
mainly in a change to the model-dependent virial radius, 
whereas for a fixed radius the dependence is significantly 
smaller. 

4. The REFLEX sample and the values of 
important model parameters 

The KL analysis is performed with a subsample of 426 
clusters of the 452 REFLEX clusters. The clusters of 
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Fig. 1. Likelihood contours (l-3a levels for two degrees of freedom) for models with h — 0.64 (left) and h — 0.80 (right). The 
remaining values of the cosmological parameters and model assumptions are the same as for the reference KL solution described 
in the main text (see also Table 1). Small stars mark the parameter values with the highest likelihoods. 
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Fig. 3. Likelihood contours for models with Qbh 2 = 0.018 (left) and Q.bh 2 — 0.026 (right). The values of the other parameters 
are given in Table 1. 



the subsample are selected according to the following cri- 
teria. They are located within comoving distances r < 
1.0/i _1 Gpc (z < 0.365). The clusters have at least 10 
X-ray source counts detected in the ROSAT energy band 
0.5-2.0 keV, X-ray luminosities L x > 2.5 x 10 42 /i~ 2 ergs" 1 
and X-ray fluxes Sx > 3.0 x 10 -12 ergs -1 cm -2 in the en- 
ergy band 0.1-2.4keV. The clusters are located in an area 
of 4.24 sr in the southern hemisphere with Declination 
< 2.5 deg, excluding galactic latitudes \b\ < 20deg and 



some additional crowded fields like the Magellanic Clouds. 
More information about the sample construction can be 
found in Bohringer et al. (2001). 

Several tests with observed and simulated data suggest 
that the REFLEX sample is at least 90% complete up to 
the lower X-ray luminosity limit given above (Bohringer 
et al. 2001, Schuecker et al. 2001). In the computation 
of the X-ray luminosities a /3-model based correction is 
applied to correct the observed flux for the missing flux 
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Fig. 4. Likelihood contours for biasing models of Sheth & Tormen (1999, left) and of Kaiser (1984, right). The former model 
yields the adopted reference solution. The values of the other parameters are given in Table 1. 
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Fig. 5. Likelihood contours for mass/X-ray luminosity relations with A = —20.055, a = 1.652, a e g = 19% (left) and A — 
— 18.320, a — 1.538, cr e g = 25% (right); definition of A and a as in RB02. The values of the other parameters are given in 
Table 1. 



outside the observational apperature radius. The correc- 
tion is on average about 10% and the correction procedure 
has been tested with Monte-Carlo simulations (Bohringer 
et al. 2002). For the cosmic ^-corrections we use a re- 
fined Raymond-Smith code, where the X-ray temperatures 
Tx are estimated with the Lx-Ix relation of Markevitch 
(1998) without 'cooling flow' corrections (ignoring higher- 
order effects caused by the fact that the latter relation was 
obtained for non-total X-ray luminosities). 

For the transformation of the X-ray luminosity limit 
-kmin(z) to the corresponding mass limit defined within 
an EdS geometry (see Sect.||) the empirical mass/X-ray 
luminosity relation of RB02 is used with the particular 
parameter values A = —20.055 and a = 1.652 (following 
their notation). They measured the masses of 106 galaxy 
clusters with redshifts generally below z = 0.1 mainly with 
ROSAT PSPC pointed observations and gas temperatures 
as published mainly from ASCA observations (assuming 
spherical symmetry, isothermality and hydrostatic equilib- 
rium). The sample is large enough to give a statistically 
representative summary of the local mass/X-ray luminos- 
ity relation of massive clusters and a rough but useful es- 
timate of its intrinsic scatter (see below). In order to test 
the sensitivity of the cosmological results on the represen- 
tativity of the sample used to derive the mass/X-ray lu- 



minosity relation, we alternatively work with an empirical 
relation derived from a smaller subsample of 63 clusters of 
RB02 which is, though strictly flux-limited, significantly 
dominated by clusters with cooling flow signatures. 

Having converted the X-ray luminosity limit in the way 
described abo ve, the resulting mass limit is transformed 
with Eq. (A. 7) in Appendix^ into a mass system which is 
consistent with the actual values of the cosmological pa- 
rameters. This transformation is necessary because the av- 
erage density contrast used to define the maximum radius 
for the mass integration depends on the values of the cos- 
mological parameters (see, e.g., Lahav et al. 1991). After 
a further transformation using (Eq. |A.8| ) these masses are 
in the system introduced by Jenkins et al. (2001) so that 
Eq. (||) can be integrated to give the expected average 
number of clusters. 

To be more specific, the computations of the two trans- 
formations are based on a general relation between cluster 
mass and density threshold (see Eq. |A.3| ). We further as- 
sume that the redshift-dependency of the critical density 
contrast might introduce a possible redshift-dependency 
of the cluster masses virialized at z (see, e.g., Eq. A. 7) 
which is of the form M(z) ~ M(0)/E(z) with the E{z) 
function already introduced in Sect. || (see also Mohr et al. 
2000). Finally, the Jenkins et al. mass function is defined 
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relative to the O m -dependent average mass density. The 
transformation of the masses (still defined relative to the 
critical - though cosmology-dependent - average density 
contrast) to the masses used for the Jenkins et al. mass 
function is accomplished via (A.8). For a given mass den- 



sity profile (see below) all these mass transformations can 
be computed in a straightforward way, and they have to be 
re-computed every time when a new cosmological model 
is tested in the likelihood analysis. 

The integration of the mass function includes a con- 
volution which takes into account (i) the intrinsic scatter 
of the mass/X-ray luminosity relation (possibly caused by 
'cooling flows' and cluster mergers) estimated to be about 
a M = 20% in mass (see below), and (ii) the random flux 
(luminosity) errors of a L = 10-20% of the REFLEX clus- 
ters (Bohringer et al. 2001) where the different error com- 
ponents are assumed to add quadratically. The intrinsic 
scatter is computed with the observed total mass scatter 
of about 50% of the empirical mass/X-ray luminosity rela- 
tion and the individual mass measurement errors of about 
30% (times a factor 1.5 described below) obtained with 
the total sample of 106 clusters of RB02. 

However, the mass measurement error of about 30% is 
formal and assumes, e.g., a constant cluster temperature 
profile and a symmetric mass distribution, which is likely 
not the case. For example, Evrard et al. (1996) studied 
the accuracy of X-ray mass estimates using gasdynamic 
simulations and found for f2 m = 1 and a critical density 
contrast of 250 a ratio of 1.15 for the isothermal to the 
non-isothermal mass estimates and a related la random 
mass error of 36%. Moreover, the Chandra observation of 
the elliptical cluster RBS797 (Jetzer et al. 2002) give dif- 
ferences between the masses derived under the assumption 
of spherical and ellipsoidal cluster shapes ranging from 
10% (oblate) to 17% (prolate). More realistic mass errors 
could thus easily be larger than the given formal error 
by a factor of 1.5. Taking this factor into account we get 
the above-mentioned maximum estimate of the intrinsic 
scatter of <7m = 20%. 

The effective scatter in mass, a c g = \J o~\i + {<tl/o) 2 , 
includes also the scatter introduced by the flux errors of 
the REFLEX clusters. The range of expected <tl values 
from 10 and 20% yields, for gm = 20%, values between 
(Jeff — 21% and 23%. In order to estimate the effects re- 
lated to our incomplete knowledge of om contributing to 
(Jeff on the cosmological constraints, we finally assumed 
in our tests cr ff values ranging from 19% to 28% with a 
default value of 25%. 

For all mass transformations the Navarro et al. (1997) 
mass density profile is used with a redshift and mass- 
independent concentration parameter of c = 5 typical for 
X-ray clusters of galaxies (see simulations of Navarro et 
al. 1997 and, e.g., the observed average value c = 5.2 ob- 
tained by Allen et al. 2002). The sensitivity of this specific 
choice on the final results is tested by using alternatively 
c = 4 and c = 6. We thus assume that the REFLEX clus- 
ters do not show any significant evolution up to z = 0.3 
as suggested by the redshift-independent distribution of 



the comoving number densities of the REFLEX clusters 
(Paper I, see also Gioia et al. 2001 and Rosati et al. 2002). 
Therefore, all mass transformations derived in Appendix^ 
and model mass functions are evaluated at the formal red- 
shift z — 0.0 where the empirical mass/X-ray luminosity 
relation was referred to, or alternatively z = 0.05, that is 
the mean redshift of the cluster sample used to determine 
the relation. 

The REFLEX clusters are counted in 360 cells (spheri- 
cal coordinates) : 6 angular bins in Right Ascension, 6 bins 
in Declination, and 10 bins along the comoving radial axis. 
With standard linear algebra codes (Press et al. 1989) we 
compute the eigenvectors and eigenvalues of the whitened 
correlation matrix R' . Tests of the KL method based on 27 
independent REFLEX-like cluster samples extracted from 
the Hubble Volume Simulation ensure the correct handling 
of the data (see Appendix of Paper II) . The higher order 
KL modes obtained for the present investigation show an 
increasing number of zero-crossings and a decreasing am- 
plitude (by definition). The same behaviour is found in 
our previous KL analysis (see Figs. 1 and 2 in Paper II) 
and is thus not illustrated again. 

Finally, it should be noted that the model covariance 
matrix C is not diagonal unless the fiducial model used to 
compute the KL eigenvectors is identical to the model used 
to compute R' . However, we found that a change to an- 
other cosmology (e.g., EdS) has a small effect on the eigen- 
vectors, slightly broadening the likelihood contours, but 
leaving the location of the likelihood maximum almost un- 
changed, as first noted in Vogeley & Szalay (1996). Tests 
show that for different initial values of fi m and as, the 
KL solutions converge to similar Q m and as values which 
we then adopted as our fiducial cosmology (see Sect.|I|). 
Hence, the fiducial cosmology is close to our final result 
(Sect. [3j) so that the effects mentioned above are negligible. 

5. Results 

Table 1 summarizes the basic results for f2 m and as includ- 
ing their formal la standard deviations obtained with the 
REFLEX cluster sample for different priors under the gen- 
eral assumption of a non-evolving cluster population and 
geometrically flat cosmologies. Typical likelihood contours 
are shown in Figs, [j] to |[ They illustrate the sensitivity of 
the analysis on important cosmological parameters. For 
the comparison of different KL solutions we use as refer- 
ence the results obtained with h = 0.70, O^/i 2 = 0.020, 
rtg = 1.0, z — 0.0, a c g — 25%, an empirical mass/X-ray 
luminosity relation with the parameters A = —20.055, 
a = 1.652 (denoted by M106), and the biasing model of 
Sheth & Tormen (1999), yielding (Fig.| left and the 8th 
row in Table 1) 



Q - 341 + 031 

"m — o.oii -0.028 ' 



rr — n 71 1 +0-039 

as - U-'-li_o.o3i 



(4) 



The la random errors do not include cosmic variance. 
Compared to the results obtained by utilizing the cluster 
fluctuations only (Paper II) , the present random errors are 
a factor of 2.0 smaller for £l m and a factor of 8.6 smaller for 
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( B „- <B „ > )/^ <B ? > 

Fig. 6. Histogram of the KL coefficients of the REFLEX clus- 
ter sample obtained with the reference solution (Eq.^), and 
superposed unit-variance Gaussian profile shifted by — 0.14er in 
order to get a better representation of the observed data. The 
coefficients are normalized by the eigenvalues of the REFLEX 
sample correlation matrix. 



er 8 . This illustrates the importance of including the clus- 
ter abundance as a second criterion because it appears to 
be the main driver of the high accuracies reached in the 
present investigation. However, the results (Q) are more 
sensitive to systematic errors compared to those solely ob- 
tained with the fluctuation analysis. Therefore, we have to 
evaluate the systematic errors introduced by the specific 
values of the priors and model assumptions in more detail. 

As expected, small systematic differences (Af2 m , Acr 8 ) 
are found when the Hubble constant h (Freedman et al. 

2001) and the spectral index of the initial scalar fluctu- 
ations ns (see references given in Sect.Q) are changed 
within their respective la ranges around the priors of the 
reference solution (see Figs.jl], 0). They are of the order of 
the la random errors given in (|4|). 

The systematics introduced by the fi^/i 2 ranges of, 
e.g., Pryke et al. (2002) appear comparatively small (see 
Fig.||). A larger effect is expected when the full range 
0.0095 < D, b h 2 < 0.023 obtained with standard Big Bang 
nucleosynthesis would have been used. However, the mea- 
surements of Deuterium in quasar absorption systems if 
primordial gives fif,/i 2 = 0.020±0.0015, which is quite con- 
sistent with recent CMB estimates (see review of Sarkar 

2002) and thus provides a good argument for the smaller 
confidence range adopted in the present investigation. 

In order to test the stability of the results with re- 
spect to the biasing model two extreme cases are selected 
(see Fig.§. The model of Sheth & Tormcn (ST, 1999) al- 
though motivated by a Press-Schechter like prescription is 
calibrated with dark matter simulations. The pure statis- 
tical biasing model of Kaiser (KA, 1984) is mainly based 



on the Gaussianity of the cosmic matter field. On the large 
scales studied in the present investigation we do not expect 
large differences between the two models, as supported by 
the results obtained with the KL analysis (see Fig. [| and 
Table 1). 

More technically, we measured the systematics intro- 
duced by the assumption that the empirical mass/X-ray 
luminosity relation of RB02 is determined alternatively at 
the formal sample mean z — 0.05 and not at z = 0.0 
yielding Afl m = -0.004 and Aa& = +0.007 (no fig- 
ure). We further tested the sensitivity of the cosmologi- 
cal parameters on the effective scatter of the mass/X-ray 
luminosity relation by alternatively assuming instead of 
cr c ff = 25 % the effective scatter of a c s — 19% and 28 % 
yielding Afi m = -0.021, Act 8 = +0.087 (compare Fig.| 
left (25%) with Fig. | left (19%)) and Afi m = +0.079, 
Acr 8 = —0.155, respectively (no figure). We also tested 
the sensitivity of the cosmological parameters on the cho- 
sen mass/X-ray luminosity relation using instead of the 
default relation obtained with the extended RB02 sample 
of 106 clusters (denoted M106: A = -20.055, a = 1.652) 
a relation obtained with a strict flux-limited sample of 63 
clusters (denoted M63: A = -18.320, a = 1.538) yield- 
ing Afi m = -0.060 and Aa 8 = +0.056 (compare Fig.| 
left (M106) with Fig.| right (M63)). Finally, we tested 
the effect of a value of the concentration parameter of the 
NFW mass density profile different from the default value 
c — 5. For c = 4 and c = 6 we obtained Afi m = +0.007, 
Aa$ = +0.014 and An m = -0.002, Aer 8 = -0.012, re- 
spectively (no figure). 

For the correct determination of the confidence ranges 
caused by the combined effect of all uncertainties of the 
priors and model assumptions one has to analyse the 
complete 8-dimensional parameter space. The KL method 
is, however, quite computer-intensive so that the errors 
could only be evaluated in the following simplified man- 
ner. The approximate combined effect (conservative up- 
per limit) is estimated by summing up the individual sys- 
tematic errors. We obtain the maximum systematic errors 
of Ail m =ta.lf 7 and Act 8 =±0.243 ■ A correlation of the 
systematics evaluated above is, however, quite unlikely. 
More realistic error estimates assume that the systemat- 
ics are uncorrelated and that the squared A's can be re- 
garded as the variances of error distributions assumed to 
be Gaussian. In this case the errors add quadratically and 
we obtain the final systematic errors 



an 



+0.087 
-0.071 i 



+0.120 
-0.162 • 



(5) 



Finally, we tested whether the Gaussianity assumed 
throughout the likelihood analysis is supported by the 
present measurements. Although the histogram of the nor- 
malized KL coefficients (Fig.|(j) shows a small systematic 
shift of 0.14cr and some excess at large positive fluctuations 
relative to a zero-mean, unit-variance Gaussian profile, the 
latter distribution and the un-shifted observed distribu- 
tion must be regarded as selected from the same parent 
distribution when tested on the 3a confidence level (KS- 
test). The present measurements thus support the funda- 
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Table 1. Constraints on Q m and as and their la errors (without cosmic variance) obtained with the KL method for the REFLEX 
cluster sample assuming a flat space, h Hubble constant. Qth 2 : baryon density, ns spectral index of initial scalar fluctuations, b 
biasing model: ST Sheth & Tormen (1999); KA Kaiser (1984). a c g relative effective scatter of mass of the empirical mass/X-ray 
luminosity relation given in percent. M/L empirical mass luminosity relation: M/L = 106 corresponding to A = —20.055, 
q = 1.652 (obtained in RB02 with 106 clusters); M/L = 63 corresponding to A — —18.320, a — 1.538 (obtained with 63 
clusters). Q m , as and their ±1<t errors (no cosmic variance); c concentration parameter of the NFW profile. Af2 m and Aas'- 
systematic differences in Q m and as between actual KL solution and reference solution. 
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Fig. 7. Likelihood contours obtained with the REFLEX clus- 
ter abundances only (center at Sl m = 0.32 and as = 0.71). The 
comparison with Fig.^j (left) obtained by using both cluster 
abundance and fluctuations shows that the KL method breaks 
the degeneracy between erg and fi m . Note the different scales 
of the diagrams. 



mental assumption of Gaussianity, although less signifi- 
cant compared to our first findings (see Paper II). 

6. Discussion and conclusions 

In the present investigation cluster number counts and 
their spatial fluctuations are analysed for the first time 
simultaneously in order to constrain cosmological models. 
The test utilizes the complementarity of clustering and 
abundance of galaxy clusters and allows us to fully exploit 
the cosmological potential of the REFLEX survey of X-ray 
clusters (but see below). 



For spatially flat cosmologies as suggested by all recent 
CMB measurements (Sect. |l|) and for a non-evolving clus- 
ter population as suggested by the z-independent comov- 
ing REFLEX cluster number densities (Paper I), we obtain 
Q m = 0.341 for the cosmic matter density and as = 0.711 
for the normalization of the matter power spectrum. The 
random errors range between 8-9% for Q m and 4-5% for 
erg, not including cosmic variance. Systematic errors can 
be identified and quantified so that the whole process ap- 
pears to be well-controlled. Under the assumption that the 
systematics add quadratically the total systematic errors 
are found to be about 2.5 and four times larger than the 
random errors of O m and as, respectively. 

Compared to the results obtained with the fluctuations 
only (Paper II), the combination of abundance and fluc- 
tuation measurements reduces the random errors of O m 
by a factor of about two and of as by a factor of about 
nine. The main driver of the accuracy of the present re- 
sults thus appears to be the cluster abundance, especially 
for as (but see below). 

In contrast to many other cosmological parameter es- 
timations we tested the assumed functional form of the 
sample likelihood function. As in Paper II we could verify 
Gaussianity on large scales - a fundamental property of 
the cosmic matter field. In this respect the present analysis 
appears to be fully consistent and self-contained. 

The KL likelihood contours (Figs. [I] to |||) show some 
correlations between fi m and as , probably indicating that 
the primordial part of the power spectrum is still less con- 
strained by the REFLEX data. The correlations are, how- 
ever, restricted to comparatively small (Sl m ,as) ranges so 
that it does not make much sense to determine an VL m -as 
relation from the data. Therefore, the degeneracy between 
Q m and as can be regarded as broken. The breaking of the 
degeneracy by the KL method using both cluster abun- 
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Fig. 8. Comparison of the reference KL solution (f2 t ot = 1, O m = 0.341, Q b h 2 = 0.020, h = 0.70, n s = 1.0, a s = 0.711, 
(j ff = 25%, ST biasing) with REFLEX observations. Left: Power spectrum obtained with the flux-limited REFLEX sample 
(Paper I, points with la error bars including cosmic variance) and the KL solution (continuous line). The theoretical model 
takes into account the effects of the different volumes covered by the present KL and the former Fourier analysis, includes the 
effects of the baryons, and is transformed into redshift space using the nonlinear model described in Paper I. In order to make 
the measured power spectrum less crowded adjacent power spectral densities and their errors were averaged. Right: Redshift 
histogram of the REFLEX subsample used here (steps) and the KL solution (continuous line). 



dance and fluctuations is nicely seen when the KL refer- 
ence results (Fig. |[ left) are compared with the likelihood 
contours obtained with the REFLEX cluster abundance 
only (Fig. |). 

Generally, likelihood analyses which optimize simul- 
taneously two or more almost independent observational 
quantities like the power spectrum and the cluster abun- 
dance, try to find the best compromise between the in- 
dividual observables and not necessarily the best fit of 
the different components. However, Fig. || shows that the 
best KL solution obtained with the REFLEX data is con- 
sistent with both the redshift histogram and the redshift- 
space power spectrum obtained with a standard Fourier 
method (Paper I). Note that the power spectrum was mea- 
sured within a comparatively small volume to reduce cor- 
relations between different fluctuation modes. A further 
argument supporting the reliability of the present cosmo- 
logical constraints is given by numerical fits of the am- 
plitude of the observed power spectrum of three volume- 
limited REFLEX subsamples using ACDM models which 
yield as — 0.70 (Paper I). This result is at variance to 
the standard normalization of the ACDM model, but in 
very good agreement with the results of the present KL 
analysis. 

The value of the cosmic matter density obtained with 
REFLEX is in good agreement with other recent inde- 
pendent measurements (see also Turner 2002). For ex- 
ample, Pryke et al. (2002) measured with DASI a mat- 
ter density of VL m = 0.40 ± 0.15 for h = 0.72 ± 0.08 
and Sievers et al. (2002) measured with CBI Q m = 



0. 62 ± 0.22 (weak priors). Netterfield et al. (2002) ob- 
tained for flat geometries Q, m — 0.33 ± 0.05, combining 
COBE-DMR, BOOMERANG, large-scale structure, and 
supernovae type la data. Sievers et al. (2002) found from 
the combination of COBE, BOOMERANG, MAXIMA- 

1, DASI, CBI, large-scale structure surveys, Hubble pa- 
rameter determinations, and supernovae type la the value 
f2 m = 0.32 ± 0.06. With the galaxy power spectrum of 
the 2dF Galaxy Redshift Survey and a new compila- 
tion of CMB data (COBE, BOOMERANG, MAXIMA-1, 
DASI, VSA see Scott et al. 2002, CBI see Mason et al. 
2002), Percival et al. (2002) obtained for a flat cosmol- 
ogy fl m = 0.313 ± 0.055 (all la errors without cosmic 
variance). Therefore, a significant number of independent 
applications appear to converge to a low matter density. 

The linear theory as obtained with REFLEX is lower 
than the COBE normalization and most estimates from 
cluster abundance as summarized in, e.g., Table 6 in RB02. 
Note that values of the effective scatter a c s of the mass/X- 
ray luminosity relation larger than assumed here would 
increase this descrepancy even more. However, small ag 
values were also obtained by Markevitch (1998) using a 
local cluster sample with ASCA X-ray temperatures and 
ROSAT X-ray luminosities, by Borgani et al. (2001) us- 
ing the ROSAT Deep Cluster Survey, by Seljak (2001) 
from the empirical mass/X-ray temperature relation of 
nearby clusters, by Viana, Nichol & Liddle (2002) using 
SDSS/RASS data and the X-ray luminosity function of 
REFLEX clusters of galaxies, and by Bahcall et al. (2002) 
using clusters selected from SDSS data. In this context 



10 



Schuecker et al.: The REFLEX Galaxy Cluster Survey 



it is also interesting to note that the KL results are in 
very good agreement with the new constraints obtained by 
Lahav et al. (2002) from their simultaneous analysis of the 
amplitudes of the fluctuations probed by the 2dF Galaxy 
Redshift Survey and COBE, BOOMERANG, MAXIMA- 
1, DASI data, and with the constraints obtained with 
the KL analysis of the early Sloan Digital Sky Survey 
by Szalay et al. (2001) obtained under similar priors (as 
described in Lahav et al.). There is thus increasing ob- 
servational evidence for a as value lower than previously, 
although some recent weak lensing results show a tendency 
for a higher cr 8 for given f2 m (e.g., Hoekstra et al. 2002, 
van Waerbeke et al. 2002). Small a s values help (James 
S. Bullock, private communication) to reduce the appar- 
ent problem of the ACDM model with the standard nor- 
malization to overpredict the number of low-mass satellite 
galaxies (Klypin et al. 1999, Moore et al. 1999). 

In all the cosmological studies mentioned above includ- 
ing the present, various model corrections might change 
the final results by 10% to 30% so that more detailed anal- 
yses of possible sources of systematic errors are needed to 
get below this error level. Our tests for systematic errors 
show that the largest systematic errors are introduced by 
the mass/X-ray luminosity relation (as expected). In order 
to use the full potential of REFLEX-like cluster surveys 
it is thus important to have more precise measurement 
of this function, i.e., its shape, intrinsic scatter, redshift- 
dependency, etc. which are expected to be provided by de- 
tailed observations with the Chandra and XMM satellites. 
Future REFLEX papers will test carefully the residual ef- 
fects of possible evolution of nearby clusters, the cosmo- 
logical constant and related quantities (H. Bohringer et 
al., P. Schuecker et al., in preparation). 

Regardless of these interim restrictions, the present in- 
vestigation illustrates that large and fair samples of X- 
ray clusters of galaxies give quite clean measurements of 
the cosmological parameters. In addition to the above- 
mentioned quality of as constraints traditionally obtained 
with clusters of galaxies, the la random error for the mat- 
ter density of a n , n = 0.030 obtained with the REFLEX 
data is already close to the random error of an m — 0.010 
expected to be attainable with the CMB Planck Surveyor 
satellite and similar to an m = 0.015 attainable with the 
SNfactory plus SNAP supernovae satellite projects (see, 
e.g., Hannestad & Mortsell 2002). Moreover, the cosmo- 
logical constraints obtained with cluster data have de- 
generacies different from high-z supernovae and CMB 
anisotropics (Holder, Haiman & Mohr 2001). Therefore, 
galaxy clusters can play a significant role in high preci- 
sion measurements of cosmological parameters. 
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Appendix A: Mass transformations 

In this appendix we discuss mass transformations which 
are necessary when masses Ma are defined within radii r& 
which correspond to different average density contrasts A 
and when average contrasts are related to different back- 
ground densities. We are particularily interested in the 



relation between the masses used in the empirical mass/X- 
ray luminosity relation of Reiprich k, Bohringer (2002) and 
the masses used in the theoretical mass function of Jenkins 
et al. (2001). To motivate the mass transformations we 
first assume an isothermal sphere at redshift z = and 
define a virial radius by 



?~vir = 



Afvm 



(47r/3)A cj0c 



1/3 



(A.l) 



with the average density contrast A c which is 187r 2 for the 
EdS case, and the critical EdS density p c = 3Hq/8ttG. 
The virial relation then reads GAfyiR = Q-TxryiR, with 
Tx the temperature of the X-ray emitting gas and a 
a structural constant. If massive clusters are simply re- 
scaled versions of low mass clusters (e.g., Neumann & 
Arnaud 2001) then a is a universal constant and one ob- 
tains 



VIR 



/a\3/2 1 



1 



1 



3/2 



(A.2) 



The assumption of an isothermal sphere mass distribu- 
tion thus yields for a given redshift and temperature the 
mass-density contrast relation A/vir/Ma = (A/A c ) 5 . 
Deviations from isothermality changes this relation as can 
be seen in Horner et al. (1999) who found from spatially 
resolved temperature measurements the empirical rela- 
tion Mvir/Ma = (A/A c ) - 266 , consistent with the mass 



density profile p(r) 



In order to be consistent 



with our reference to high-quality dark matter simula- 
tions we assume a Navarro, Frenk &: White (NFW, 1997) 
profile. In this case we have to deal with a mass ratio 
Mvir/A^a = n(A c , A) which can only be obtained nu- 
merically using 



Af Al 
M A2 



Ai 
A 2 



rAi 

fA 2 



1,1 , ca, /r 2 oo 

M Al _ ln(l + cr Al /r 200 ) - 1+cr \ JrwQ 



A/ A2 



ln(l + crA 2 /r 2 oo) 



cr &2 /r 2 oo 
l+crA 2 /r 20 



(A.3) 



(A.4) 



where c is the concentration parameter which we assume 
to be independent of mass and redshift. Iterative solu- 
tion for r^ ± /r^ 2 yields the mass ratio n where the nor- 
malization radius r2oo of the NFW profile cancels out. 
Due to its definition via a mass ratio, II has the property 
n(Ai,A 2 ) -n(A 2 ,A 3 ) =n(Ax,A 3 ), yielding 



Afvm(0) 
M 2 E o d o S (0) 



n(A c (0),A^ dS (0) = 18tt 



(A.5) 



The arguments indicate that the quantities are deter- 
mined at redshift z — 0. Note that the density contrast 
A^ ds (0) approximately defines the same mass Af^f as 
used in the empirical mass-luminosity relation of Reiprich 
& Bohringer. 

With the function E(z) defined in Sect.|^, a possi- 
ble redshift-dependence might be introduced by p c (z) ~ 
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E 2 (z) and A c (z) = 18tt 2 +82 p m {z) - l]-39 [Q. m {z) - if 
with n m (z) = ^m(l + z) 3 /E 2 (z) for spatially flat geome- 
tries (Bryan & Norman 1998). For a given temperature we 
thus obtain 



so that the relation between a virialized mass at redshift 
z and the mass as defined in the empirical mass/X-ray 
luminosity relation becomes 



Mvm(s) = M 2 E d s (0) • n(A c (z),Af dS (0)) • -L . (A.7) 



Finally, we have to relate Mvir(z) to the mass as defined 
in the universal mass function of Jenkins et al. (2001). 
Here, the masses are defined by the radius of the spheri- 
cal overdensity of 187r 2 with respect to the fi m -dependcnt 
background density. Denoting this mass by M SO (i8o)( z ) 
we have 



where f = r SO (i80)/ r viR is obtained iteratively from the 
NFW mass density profile using 



Mvir(z) = Mvir(0) 



n(A c (z),A e (0)) 
E(z) 



(A.6) 



M SO (i80){z) = M vm (z) Q m (z) r 3 , 



(A.8) 



ln(l + cf)-l + 



1 



= n m (z) ln(l + c)- 



c 



.(A.9) 



1 + cf 



1 + c 



In the main text a non-evolving cluster population is as- 
sumed so that the relevant equations are evaluated either 
at z = 0.0 or alternatively at z = 0.05. 



